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Abstract. Numerical simulations of kinematic dynamo action in steady and three-dimensional ABC flows are 
presented with special focus on the difference in growth rates between cases with single and multiple periods of 
the prescribed velocity field. It is found that the difference in growth rate (apart from a trivial factor stemming 
from a scaling of the rate of strain with the wavenumber of the velocity held) is due to differences in the recycling 
of the weakest part of the magnetic held. The single wavelength classical ABC-flow experiments impose stronger 
symmetry requirements, which results in a suppression of the growth rate. The experiments with larger wave 
number achieve growth rates that are more compatible with the turn-over time scale, by breaking the symmetry 
of the resulting dynamo generated magnetic field. Differences in the topology in cases with and without stagnation 
points in the imposed velocity field are also investigated, and it is found that the cigar-like structures that develop 
in the classical A=B=C dynamos, are replaced by ribbon structures in cases where the flow is without stagnation 
points. 
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1. Introduction 

Dynamo theory aims at explaining how magnetic fields 
can be generated and sustained by fluid motions, in the 
presence of dissipation. Studies of dynamo action in astro- 
physical objects are usually connected with the question 
of th e existence of fast dynamos (Vainshtein & Zeldovich 
|l972 ) . The fast dynamo problem refers to the generation of 
magnetic fields in the limit of infinite magnetic Reynolds 
number Re m , and has been extensively reviewed by e.g. 
Childress & Gilbert ( |1995| ). 

In the simplest formulation of the kinematic dynamo 
problem, an initially weak magnetic field is amplified due 
to the motions of a prescribed, stationary velocity field. 
Flows that have chaotic transport properties are prime 
candidates for fast dynamos. They have the ability to am- 
plify the magnetic field exponentially by stretching the 
magnetic field lines. One such class of flows that is often 
invoked in fast dynamos studies is the ABC flows (named 
after Arnold, Beltrami & Childress). A first detailed re- 
port of the kinematic dyna mo act ion of ABC flows may 
be found in Dombre et al. ([1986D , who examined espe- 
cially the internal symmetries of the velocity field. There 
is, however, no guarantee that such symmetries are pre- 
served in the eigenmodes of the magnetic field. Some dis- 
cussion of possible symmetry breaking in the eigenmodes 
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ma y be fo und in Galloway et al. ( 1984 ) and in Galanti et 
al. (1993). They found no evidence of symmetry breaking 
when the magnetic Reynolds number is limited to very low 
values, although symmetries are broken for higher values 
of the magnetic Reynolds number. The presence of sym- 
metry breaking was also shown by Galloway et al. ( 1984 ) 
when the wavenumber of the periodic ABC flow is equal 
to 3 and 4. The conclusion was that the physical reason 
for the symmetry breaking is unclear and requires further 
investigation. 

Dynamo action in ABC-flows has also been studied to 
determine the growth rate of the magnetic field for differ- 
ent values of the magnetic Reynolds number (Arnold et 
al. |!983t Galloway et al. |1984| Lau & Finn |1993| ). Results 
for the case with unit wavenumber, for magnetic Reynolds 
number up to Re m = 1600 (Dorch 2000| ) showed that the 
growth rate saturates at a value close to 0.07. Galanti 
et al. (1993) found that doubling the wavenumber of the 
flow leads to a sharp increase of the growth rate, while for 
larger wavenumbers the growth rate scales with the strain 
rate of the flow (i.e., in proportion to the wave number). 

In the current paper we address some of the issues 
mentioned above by means of numerical kinematic MHD 
simulations. The dynamics of the magnetic structures in 
physical space is examined in connection with the am- 
plification mechanism which drives the dynamo action in 
ABC flows. We also investigate the physical reason that 
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causes the increase of the growth rate of an ABC flow for 
wavenumbers larger than unity. To achieve this we follow 
the motion of magnetic field lines in time and demonstrate 
that the transportation of the weak part of the field is a 
key factor of the dynamo process. 

The numerical aspects of the simulations are described 
in the following section. The growth rates produced by 
fast kinematic dynamo action of ABC flows are presented 
in Section |j. Visualizations of the magnetic structures 
and field lines in physical space are shown in Section 
^. Dynamo action produced by flows without stagnation 
points is studied in Section ||. Section [| contains the overall 
conclusions of the numerical simulations presented here. 

2. The simulation 

Our aim is to study kinematic dynamo action in ABC 
flows. The time evolution of B is governed by the (dimen- 
sionless) induction equation: 



f)Ti 1 

— =Vx(uxB) + — V 2 B, 
at Re™ 



V B = Vu = 0, 



(1) 
(2) 



The velocity field is chosen to have the form of an ABC 
flow: 

u abc = ^4(0, sinfca;, cosfcx) + i?(cos ky, 0, s'mky) (3) 
+ C(sin kz, cos kz, 0) , 

where k is the wavenumber of the flow. The special case 
with parameters A=B=C=fc=l is referred to as the "nor- 
mal" ABC flow dynamo. In the normal case there are 
two windows of dynamo action (where the growth rates 
are positive): A closed one in the interval Re m =8.9-17.5, 
and an open one for magnetic Reynolds numbers above 27 



(Galloway et al. |1984|). 

Eq. (^) is solved numerically on a staggered mesh us- 
ing array valued functions for evaluating the space and 
time deri vativ es. A third-order predictor-corrector method 
(Hyman 1979 ) is used for the time-stepping. The applied 
numerical scheme (by Galsgaard, Nordlund and o thers , 
see Galsgaard & Nordlund [l997| Nordlund et al. |l992|) 
ensures Eq. (||) . Periodicity of 2ir is assumed over a three- 
dimensional computational domain. The initial magnetic 
field is a random perturbation with an amplitude of 10~ 5 . 

3. Growth rates 



Previous numerical simulations (e.g. Arnold et al. |1983 



Galloway et al. 1984 



1993; Dorch 2000) 



Galanti et al 

have shown that the normal (A=B=C=fc=l) ABC flow 
very likely acts as a fast dynamo with a growth rate close 
to 0.077 (corresponding to a time-scale of ss 0.25 turn- 
over times) . When the wavenumber k is larger than unity 
(where k = 1 corresponds to the largest scale of 2tt in the 
computational domain) the growth rate scales approxi- 
mately with the rate of strain. These results are confirmed 



in our simulations and we continue the study of the growth 
rate of the magnetic field for higher values of Re m . 

We first discuss the case with A=B=C=1 and k = 2 
before turning our attention to higher wavenumbers: As 
shown in Fig. (]]]) the magnetic field is amplified exponen- 
tially without oscillations. Initially, several modes may be 
present, associated with oscillatory growth of the mag- 
netic field (as in the case of a normal ABC flow, Dorch 
|2000 ), but eventually the mode with the highest growth 
rate is the only one that remains in the solution (i.e. the 
dominant mode). 




20 30 40 
Time 

Fig. 1. Temporal evolution of the magnetic energy for the 
kinematic ABC dynamo with A=B=C=1 and k = 2 for Re m 
liifi and Re m = 1600. 



For low Re m the growth rate increases with Re m and 
reaches its maximum value at Re m = 12. This value of 
the magnetic Reynolds number belongs to the first dy- 
namo window for the normal ABC flow dynamo, and cor- 
responds to a maximum growth rate of the magnetic field 
in this window (Galanti et al. 1993). With k — 2 the lowest 
Re n 



possible for dynamo action is found to be approxi- 
mately equal to 2. For large Re m the growth rate decreases 
sli ghtly (in contrast to the low Re m results of Galanti et 
al. 1993| ) but always remains positive and bounded away 
from zero (Fig. [|). The largest value of Re m considered 
here is Re m = 1600. In the high Re m regime the growth 
rate saturates at a value close to 0.3. 

The evidence indicates that the ABC flow with k = 2 
is a fast dynamo, since the growth rate seems to saturate 
at a non-zero value for high Re m . The reason for the en- 
hancement of the growth rate when k — 2 is discussed in 
the following section. 

The results for the growth rate of the magnetic energy 
at different values of magnetic Reynolds number when k = 
3 and k — 4 are shown in (Fig. 0). A comparison of these 
results with the k = 2 case (FigH3) shows that the growth 
rate of the magnetic field increases with increasing wave 
number, in approximate proportionality to k for k > 2, 
scaling approximately as 0.15 x k, while the growth rate 
is unproportionally small (ss 0.07) for k = 1. 
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Fig. 3. Two panels showing the growth rate of the magnetic energy normalized by the rate of strain for the kinematic ABC 
dynamo with k = 3 (left panel) and k — 4 (right panel). 
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Fig. 2. Growth rate of the magnetic energy (normalized by 
the rate of strain) versus the magnetic Reynolds number for 
the kinematic ABC dynamo with A — B — C = 1 and k = 2. 



4. Magnetic structures 

An experiment with Re m = 100 and a numerical resolu- 
tion of 64 3 was performed in order to study the dynam- 
ics of the magnetic structures in physical space when the 
wavenumber k of the flow is equal to 2. The computational 
domain then contains eight periodic cells, where each one 
(initially) reflects the geometrical properties of a standard 
ABC flow. The magnetic field evolves according to the 
induction equation Eq. (Q) and flux "cigars" are formed 
at the so-called a-type stagnation points (Dombre et al. 



1986). These are points (stable manifolds of the flow) with 



converging streamlines in the plane of the stagnation point 
and the field is rapidly advected to these regions (Fig. ||). 
The flux cigars are aligned along the axis of divergence 
through the a-type stagnation points. 

At the /3-type stagnation points (unstable manifolds 
of the flow) the field is weaker, and takes the form of flux 
sheets on both sides of the plane of divergence. The flux 
sheets are formed because of the converging flow along 
the axis of divergence from the a-type stagnation points. 




Fig. 4. Magnetic field strength isosurfaces (dark) and stag- 
nation points of the flow (light). The magnetic field is con- 
centrated along flux cigars centered at the a-type stagnation 
points of the flow. The stagnation points with no flux cigars 
centered at them are the /3-type points. In each cell three cigars 
form a triangle around a /3-type stagnation point and another 
one is pointing towards it. The magnetic field strength of the 
isosurfaces is 70% of the peak value in the snapshot. 



The flux cigars thus appear to be connected by flux sheets 
through hctcroclinic stream lines. 

When the symmetry breaking appears, and the fastest 
growing eigenmode takes over, the topology of the mag- 
netic field changes. The fastest growing mode does not 
follow the symmetric configuration of a standard ABC 
flow — the strong magnetic field is still elongated into flux 
"cigars" but their size and field strength vary from cell to 
cell (Fig. |). 

The rate of growth (and decay) of individual flux cigars 
is different, rather than being uniform as in the normal 
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Fig. 5. A snapshot showing the structure of the magnetic field 
when the fastest growing mode appears and the symmetry is 
not conserved. Each cell consists of flux cigars which have dif- 
ferent size and different field strength. 



k = 1 case. The magnetic field for the latter case is re- 
cycled inside its cell because of the exact symmetry. The 
time evolution of the magnetic energy follows a periodic 
(oscillatory) behavior which is associated with the con- 
servation of symmetry in the fastest growing eigenmode 
(Dorch 2000) in contrast to the dominant mode for the 
k — 2 case which is a symmetry breaking mode (Galloway 
et al. 



1984). Visualization of magnetic field lines shows 



that they are not "trapped" in each individual cell but 
there is always a "loose" weak field which is transported 
from cell to cell and is stretched up against the local flux 
cigars that are passed on the way (Fig. |^). In general, 
the increase of the magnetic energy is associated with the 
stretching, twisting and folding of the magnetic field lines 
against the local flux cigars in the /3-type planes. The field 
lines coming from a flux cigar are connected through the 
weak field to the cigars that form the triangles around the 
/3-type stagnation points. The field lines at the corners of 
the triangle are stretched and added to other flux cigars of 
neighboring cells through the diverging stream lines. The 
flux sheets that are formed on both sides of the /3-type 
plane are also folded constructively. The /3-type plane is 
actually a plane of discontinuity where the magnetic field 
lines change direction and are released to move back into 
the general flow (Fig. fy. 

5. Flows without stagnation points 

In view of the role played by the stagnation points in the 
normal ABC flow dynamo, it is natural to ask whether 
they are a necessary ingredient for fast dynamo action. 
When the parameters A, B and C are varied it is possi- 
ble for the stagnation points to be created or destroyed 




Fig. 6. Visualization of magnetic field lines and flux cigars 
when the dominant mode starts to take over. The stretching 
of the field lines at the corner of a triangle at the bottom-left 
cell is obvious. These field lines get pressed up against another 
cigar in a neighboring cell (up-right). 




Fig. 7. Folding of week flux sheets at the plane of a /3-type 
stagnation point. The field lines (top of the figure) are coming 
from a neighboring flux cigar and change direction across the 
plane. 



(Dombre et at. |l986|) . When A=5, B=2, C=2 or A=5, 
B=2, C=l (abbreviated as the 5:2:2 and 5:2:1 case, 
respectively) the flow has no stagnation points. Numerical 
simulations were performed for both cases and the results 
obtained are similar. Therefore, only the former case is 
discussed in detail. 

The 5:2:2 flow has four points where the velocity is 
close to zero (but not zero, hence there are no stagnation 
points). The disappearance of the stagnation points means 
that there are no longer any heteroclinic stream lines that 
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connect the points of minimum velocity. However, there 
are stream lines that converge to these points and stream 
lines that diverge from these points (Fig. ||). The stream 
lines of the 5:2:2 flow do not follow any specific sym- 
metry as in the normal ABC flow (which has a three-fold 
symmetry of stream lines around the stagnation points) 
and the topology of the flow becomes more complex. 




Fig. 8. The bead-shaped isosurfaces show the position of the 
points with minimum velocity of the 5:2:2 flow. Visualization 
of stream lines of the flow shows that they diverge when they 
reach the neighborhood of other points with minimum velocity. 

The importance of these minimum velocity points is 
that they coincide with local extrema of the "perpendic- 
ular convergence" of the local motion to a magnetic field 
line 



D ± ]n.B/Dt= -V 



(4) 



which, apart from diffusive effects, gives the exponen- 
tial rate of increase of the local magnetic field amplitude 
(Nordlund et al. 1994). In general, it is the gradients of 



the flow velocity that enable the stretching of field lines 
that is so essential to the dynamo mechanism, and conse- 
quently these points are points where the stretching is at 
maximum. 

5.1. Growth rate 

Numerical experiments were performed with Re m up to 
3200 and numerical resolution up to 240 3 to test whether 
the 5:2:2 case has vigorous dynamo properties in the 
kinematic regime. The initial seed field was a random per- 
turbation with an amplitude of 10 -5 . It was found that the 
magnetic field is amplified exponentially and the growth 
rate increases with Re m as illustrated in Fig. (||). 

The 5:2:1 case was also examined at high magnetic 
Reynolds number. The highest value of Re m considered 
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Fig. 9. Growth rate of the magnetic energy versus the mag- 
netic Reynolds number for the 5:2:2 flow. 



here is Re m = 3200. The growth rate seems to saturate at 
a value of about 0.67 as illustrated in Fig. (Kw. 
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Fig. 10. Growth rate of the magnetic energy versus the mag- 
netic Reynolds number for the 5:2:1 flow. 



5.2. Flux ribbons 

As mentioned above, the 5:2:2 flow has four points with 
minimum velocity that coincide with local maxima of the 
function in Eq. 

When the induction equation is evolved from a weak 
random field, flux ribbons arise close to the minimum ve- 
locity points. The thickness of the ribbons decreases as 
Re m increases, presumably scaling as Re m ' , similarly 
to the flux cigars for the standard ABC flow. The sign 
(polarity) of their magnetic field changes from ribbon to 
ribbon. Thus, field lines with opposite polarity may recon- 
nect when two ribbons are close together. The magnetic 
flux is released through the reconnected field lines (Fig. 
pd| ). However, the magnetic ribbons are not confined at 
certain points but are advected by the flow, and field lines 
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with the same polarity could also be found to approach 
each other. 




Fig. 11. The magnetic field takes the form of flux ribbons 
when A : B : C = 5 : 2 : 2. The reconnection of field lines 
is shown when ribbons with opposite polarity field are found 
close together. 




Fig. 12. A snapshot showing the merge of the ribbons (dark) 
with like polarity in the neighborhood of the minimum velocity 
points (light). 



Due to the stretching and constructive reinforcement 
of the magnetic field discussed above, the magnetic energy 
increases each time the process is repeated. 



5.3. The amplification mechanism 

The process responsible for driving the amplification 
mechanism of the mode with the largest growth is the 
constructive folding of the ribbons. 

Two of the most important properties of the 5:2:2 
flow is the exponential stretching of the field lines and the 
way in which the magnetic field is advected and the field 
lines are packed. Neighboring regions with magnetic field 
of opposite polarity are destroyed whereas regions with 
like polarity flux are conserved. Therefore the rate of field 
generation depends on the global folding of the field lines. 

The flow folds the field so as to bring together flux in 
the same direction. Thus, stretched flux ribbons with the 
same polarity come into alignment and merge. The con- 
structive folding of the sheets occurs in the vicinity of the 
minimum velocity points (Fig. |l2] ). Due to the stretching 
ability of the flow, which is maximum close to these points, 
the field lines are stretched in these regions. 

The flow advects the field so as to bring together 
stretched field lines with the same polarity. They come 
into alignment and are folded constructively. These field 
lines are continuously stretched and also advected by the 
flow, and move to the vicinity of another neighboring 
point with minimum velocity, where they fold with new 
stretched field lines of the same polarity. The folding is 
again constructive. 



6. Conclusions 

The amplification of the dynamo generated field is asso- 
ciated with a combination of stretching, twisting, folding 
and reconnection of magnetic field lines. The importance 
of the weak part of the field in the dynamo process is 
emphasized when the wavenumber of the flow increases. 
Symmetry breaking allows the transportation of stretched 
weak magnetic field from cell to cell, where it gets packed 
against the local flux structures ("cigars"). Most of the 
space in the bulk part of the flow has passive advection, 
which does not change at increasing magnetic Reynolds 
number. The evidence indicates that the invariant prop- 
erties of the flows (i.e. stretching ability) do not change in 
the limit of arbitrarily small magnetic diffusivity. 

The study of kinematic dynamo action in ABC flows is, 
from a physical point of view, simple enough since it does 
not contain the complexity of non-linear back-reaction of 
the magnetic field. A fully non-linear study of the dynamo 
problem connects the initial phase of the growth of a weak 
magnetization state with the saturation of the dynamo 
generated field. 

We have demonstrated that the fast dynamo action 
reflects the same mechanisms at work whatever the val- 
ues of the parameters (such as Re m and k) are. The main 
mechanism is the stretching of the weak magnetic field. 
This is a first step towards the understanding of the pro- 
cesses at work in a turbulent environment and during the 
saturation state of a dynamo where the flow has still a 
good grip on the weakest part of the magnetic field. 
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It is of course also of great interest to perform dynamo 
simulations in more realistic (e.g. spherical) geometries, 
and to include more physical processes such as differential 
rotation and convection. The advantage of the simulations 
presented in the current paper is that the individual mech- 
anisms are easily identified and can studied in great de- 
tail. Thus, kinematic dynamo experiments such as these 
can provide very useful clues for the general theory of as- 
trophysical dynamos. 
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